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Stable chaos in the 55Cnc exoplanetary system? 
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ABSTRACT 

The five planets discovered around the main-sequence star 55 Cnc may represent a 
case of stable chaos. By using both the Frequency Map Analysis and MEGNO we find 
that about 15 % of the systems that can be build from the nominal orbital elements of 
the system are highly chaotic. However, in spite of the fast diffusion rate in the phase 
space, the planetary system is not destabilized over 400 Myr and close encounters 
between the planets are avoided. 
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1 INTRODUCTION 

The most crowded extrasolar planetary system discovered 
so far is that around the main-sequence star 55 Cnc(= p^ 
Cancri). Doppler shift measurements strongly suggest the 
existence of five planets orbiting the star with semimaj or 
axes ranging from 0.038 to 5.901 AU (|Fischer et al.lbOOsi ') . 
The innermost and smallest planet, e, is a Neptune-mass 
object while the outer and most massive one, planet d, 
has a minimum mass of about 4 Jupiter masses. A self- 
co nsistent dynam i cal fit to the stellar wobble data performed 
by iFischer et al.l (|200a ) gives a set of orbital elements for 
the p lanets reported in Table 1. It is claimed IjFischer et al.l 
120081 ) that this nominal system of five planets described in 
Table 1 is dynamically stable at least over a timescale of 1 
Myr. We would like to point out that planets 6 and c are 
no longer in a 3:1 me an motion resonance as in the previ- 
ous solution given by iMcArthur et al.l 1120041'). In addition, 
the eccentricity of planet c in lFischer et al.l (120081) is signifi- 
cantly smaller compared to that derived bv lMcArthur et al.l 
l|2004l ). 

We have performed a detailed exploration of the stabil- 
ity of the nominal system by applying the Frequency Map 
Ana lysis (hereina fter FMA ) method I L askar et al.l 1 19921 : 
[Laskar 1993a, b; Marzari et a l. 2006) on 400 varied systems. 
These systems are obtained by varying orbital elements of 
the nominal system. The varied systems are analysed by 
adding randomly some inclination (lower than 5 deg) to 
the Keplerian orbits to make the system more realistic. The 
masses of the planets are scaled accordingly. Among all the 
systems analysed with FMA we obtain cases with large dif- 



fusion speed in the phase space suggesting a chaotic evolu- 
tion and possible instability. However, integrating a few of 
these chaotic systems over timescales of 10^ yrs, we never ob- 
tain destabilization. We, therefore, think that these chaotic 
system a re examples for stable chaos. This term was intro- 
duced bv lMilani fc Nobilil (119921 ) to indicate the peculiar be- 
haviour of asteroids with short Lyapunov times (of the order 
of some 10"^ yrs) which however show a remarkable sta bility 
over the age of the solar system (|Milani fc Nobililll992l ). For 
the 55Cnc planetary system we find that in about 15% of 
all cases the orbital elements of the four inner planets ex- 
hibit a chaotic evolution on a timescale of a few million of 
years. This stable chaos may be due to the topology of the 
phase space: Quasi-periodic solutions in the proximity of a 
chaotic island can ac t like "quasi-barriers" for the diffusion 
(|Tsiganis et al.l 120001 ) . The chaotic evolution appears to be 
confined and the system does not destabilize at least on a 
timescale of 400 Myr. This peculiar behaviour may be due to 
the closeness of the system to the planetary 3:1 mean motion 
resonance between planets h and c and secular resonances. 

To confirm that stable chaos is not limited to the nom- 
inal case but can be retrieved also in a larger region of the 
phase space, we randomly varied the orbital elements of 
the nominal case and found a similar behaviour. For a few 
selected cases, we also applied the MEGNO method (the 
acronym of Mean Exponential Growth factor of Nearby Or- 
bits) proposed by Cincotta fc Simo (2000). The results ob- 
tained by MEGNO confirm the FMA resuhs. 
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Table 1. Orbital parameters for the self-consistent dynamical fit of the 55 Cancri planetary systems. These data are taken from Table 
4 of Fischer et al. (2008). 
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Figure 1. Histogram showing the dsitribution of the diffusion 
rate within our sample of 400 55Cnc planetary systems. 



2 FMA ANALYSIS 

We have performed the FMA analysis on 400 different sys- 
tems derived from Table 1. Random mean anomalies, node 
longitudes and inclinations lower than 5° are assigned to 
each planet. All other orbital elements were taken from Ta- 
ble 1. The frequency analysis is performed over a timescale 
of 2 X 10'' yrs. The orbital elements computed w ith the nu- 
merical integrator SYMBA (IDuncan et al.lll998l ) are Fourier 
analyzed and the values of intrinsic frequencies are obtained 
over running windows. The relative changes of these frequen- 
cies are estimated and used to compute the diffusion rate in 
the phase space. We have applied FMA to the signal satcj ^ , 
the difference between the periapse longitude of planets b 
and c, and to the more conventional signal s = he + ike- 
The usual non-singular variables he = ecc * cos{tjj) and 
ke = ecc * sin{-uj) refer to the innermost planet e where 
ecc refers to its eccentricity. The diffusion rate is computed 
by using the standard deviation a of the main frequency of 
the signal computed over the running windows. Slow diffu- 
sion rates, characterized by small values of ct, mean quasi- 
periodic systems, while largo values of a imply chaotic evo- 
lution. As in Marzari et al. (2006) we measure the diffusion 
speed by the logarithmic number Vd = ~logio{a) -'rlogia^ao) 
where uo is the smallest value of a we observed in our sam- 
ple of systems. A small value of Vd means a low dispersion of 
the frequencies in the considered time interval and, hence, a 
slow diffusion speed in the phase space. Large values of Vd in- 
dicate fast changes of the system frequencies and, therefore, 
chaos. 

In FigUwe show the distribution of the diffusion speed 
measured by Vd in our sample of 400 55Cnc planetary sys- 



tems. Values smaller than 4 suggestlong term stability ac- 
cording to our previous experience with FMA, while larger 
values indicate chaos. We concentrate on systems with a dif- 
fusion speed of about 5 or larger which represent about 15% 
of the whole sample we analysed. These systems have major 
frequencies which change on short timescales. If the fast dif- 
fusion of major frequencies induces also drastic changes of 
amplitudes, in particular of planetary eccentricities, close en- 
counters between planets may occur resulting in the ejection 
of one or more planets. We will show in the next section that 
planetary eccentricities do not increase on longer timescales 
which suggets a "stable chaos" state for the Cnc55 system. 
An additional set of FMA simulations has been per- 
formed for 400 systems where we have randomly varied the 
last signifcant digit of all the orbital elements given in Table 
1. We obtain a histogram that is very similar to that shown 
in Fig[l] confirming that the behaviour of stable chaos we 
have found for the nominal case of Table 1 is extended in 
phase space. 



3 LONG TERM EVOLUTION OF THE 
CHAOTIC SYSTEMS 

For those systems with a large diffusion speed {vd > 4) we 
performed numeric al integrations over a timescale of 400 
Myr with SYMBA (|Duncan et al.lll998l ). In Figl2]we show 
for a typical run the evolution of the eccentricity of the 3 in- 
ner planets e, h, c. The eccentricity jumps reveal the chaotic 
nature of the system. The jumps are obviously correlated 
among the three planets. Major jumps occur simultaneously. 
Also the eccentricity of planet / shows simultaneously jumps 
but more moderate. All jumps, although they appear to 
be significant, do not result in close approaches among the 
planets. The planetary system is not destabilized on this 
timescale while it is chaotic. 

Fig[2]suggests that we have the case of stable chaos. The 
eccentricity evolution is characterized by jumps. It does not 
have a random walk growth but remains limited. It appears 
that the system is bouncing between limited regions. 

A different case with similar high diffusion speed Vd is 
shown in Fig[3] The semimajor axes of the three inner plan- 
ets are shown on a shorter timescale to highlight the chaotic 
evolution. Also in this case the system remains confined in 
a stable configuration over 400 Myr. 
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Figure 2. Evolution with time of the eccentricity of the three 
inner planets of the system. The top plot shows the eccentricity 
of 55Cnc e, the medium plot that of 55Cnc b and the bottom plot 
that of 55Cnc c. 
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Figure 3. Semimajor axis of the three inner planets of the system 
over a short timescale. 



4 MEGNO ANALYSIS 

Besides the FMA method which uses the diffusion rate of in- 
trinsic frequencies to measure chaos, there is a large class of 
methods which use the divergence of nearby orbits in phase 
space as a measure. Exponential divergence means chaos 
while linear divergence means non-chaotic which guarantees 
stability of the dynamical system. Stability means here that, 
in particular, no planet is ejected out of the system. A nearly 
exponential divergence does not necessarily mean instabil- 
ity in our sense over the lifetime of the system. As outlined 
above, we would qualify such a system as stable chaotic. The 
widely used Lyapunov characteristic numbers (LCNs) yield 
the necessary information about the divergence and are used 
to measure the chaocity of the system. Since the computa- 
tion of the LCNs is very time consuming, more rapid meth- 
ods were introduced which were shown to yield a very good 
estimation for t he LCNs. One of these fast LCN estimators 
is the MEGNO ijCincotta fc Sim6ll2000l') method. It wa s ap- 
plied to several exoplanetary systems ( Bois et al.ll2003l ). We 
computed the MEGNO indicator for about 50 varied sys- 
tems. They were obtained by varying only the eccentricites 
of planets e,b,c and f while taking the other orbital elements 
from the nominal system of Table 1. Eccentricites of planet 
e are varied between 0.22 and 0.26, of planets b and c be- 
tween 0. and 0.10 and of planet / between 0. and 0.04. We 
integrated the systems over 100 000 yrs. For eccentricities 
of planets 6 and c ranging from 0. to 0.02, we obtain weak 
chaos. The estimator for the maximal Lyapunov exponent 
is not linear but far from exponential. For larger eccentrici- 
ties, the indicator has a behaviour between linear and expo- 
nential. MEGNO shows, like FMA, that the four planetary 
orbits close to the nominal system are chaotic. 



5 CONCLUSIONS 

Among the possible dynamical configurations of the 55Cnc 
planetary system there is also stable chaos. About 15% of 
the systems built f rom the nominal dynamical fit given in 
Fisc her et al.l IJ2008I ) have a high diffusion speed in the phase 
space. However, the chaotic island seems to be limited in 
extent and surrounded by quasi-periodic systems. In this 
way, the chaotic systems do not increase their eccentric- 
ity to planet crossing values and they are stable over a 
long timescale. Systems with intermediate values of diffu- 
sion speed Vd as measured by FMA (~ 4) do not show this 
behaviour. They possibly populate the outer border of the 
chaotic region and their chaotic evolution is much slower. 
Systems with smaller values of Vd are possibly stable over 
timescale of lO" yrs. Additional FMA and MEGNO com- 
putations show that this behaviour is not only peculiar for 
the nominal system given by Table 1 but it is extended over 
a wider range of orbital elements for the system. The weak 
resonances responsible for the stable chaos are present also 
for different combinations of the initial orbital elements. 
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